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Q , Table 1. Characteristics of H II regions in NGC 3256. 
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. Note. — See Figure 1 for H II region labels. The velocities are heliocentric and 

the formal error in velocity for these measurements is 7 km s~^. However the velocity 
^ , is that measured at the position on the DBS slit which lies closest to the centre of 

H ' the polygonal region used in CCD photometry; the difference between the tabulated 

velocity and the velocity at the centre of the H II region can be as large as 20 km s . 
The photometry is described in § 2.1.3; corrections to the flux density compensate 
for atmospheric extinction, interstellar medium absorption in the Milky Way, and 
filter transmission losses. The total luminosity is the Ha flux over 47r steradians at a 
distance of Vsystemic/100 h km s^^ Mpc^^. 
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Table 2. Dynamical Masses in NGC 3256. 



Radius Mass Position of Radius 

arcsec kpc 10^^ Ii^^Mq 



70 10 2.5 ± 0.5 At the edge of I band envelope 

22 3 0.8 ± 0.2 At the outer edge of HII "ring" 

13 2 0.5 lb 0.1 At the turnover peaks in optical velocity curve 



Note. — The radius used in the mass calculation is given in both angular and 
linear units. The dynamical mass calculation is described in § 3.1. 
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ABSTRACT 

CCD images and spectra of ionized hydrogen in the merging system NGC3256 were 
acquired as part of a kinematic study to investigate the formation of globular clusters 
(GC) during the interactions and mergers of disk galaxies. This paper focuses on the 
proposition by Kennicutt & Chu (1988) that giant H II regions, with an Ha luminosity 
> 1.5 X 10^" erg s~^, are birthplaces of young populous clusters (YPC's ). 
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. Although NGC 3256 has relatively few (7) giant H II complexes, compared to some 

' other interacting systems, these regions are comparable in total flux to about 85 30- 



Doradus-like H II regions (30-Dor GHR's). The bluest, massive YPC's (Zepf et al. 
1999) are located in the vicinity of observed 30-Dor GHR's, contributing to the notion 



o 

' that some fraction of 30-Dor GHR's do cradle massive YPC's, as 30 Dor harbors R136. 

^ ■ If interactions induce the formation of 30-Dor GHR's, the observed luminosities indi- 

ct ■ 

cate that almost 900 30-Dor GHR's would form in NGC 3256 throughout its merger 



epoch. In order for 30-Dor GHR's to be considered GC progenitors, this number must 
be consistent with the specific frequencies of globular clusters estimated for elliptical 
5^ I galaxies formed via mergers of spirals (Ashman & Zepf 1993). This only requires that 

about 10% of NGC 3256's 900 30-Dor CHR'sharbor YPC's which survive several Gyr 
and have masses > Mrisq- 

Subject headings: galaxies: individual (NGC 3256) — galaxies: interactions — globular 
clusters: general 



1. Introduction 
1.1. Motivation for Kinematic Observations 



Based on observations of starburst activity and the large amounts of molecular gas in inter- 
acting galaxies, Schweizer (e.g. 1988) suggested that merger events provide an ideal environment 
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for the formation of massive star clusters. Theoretical support that these would evolve into globu- 
lar clusters (GC's) was provided by Elmegreen & Efremov (1997) who argued that massive bound 
clusters preferentially form in a high pressure environment such as that created during a starburst. 
Starbursts frequently are observed in galaxy-galaxy interactions and mergers, presumably enhanced 
by tidally induced shocks (Larson 1987). Other suggestions of how interactions may provide the 
dynamical phenomenon necessary to produce young populous clusters (YPC's ; also known as super 
star clusters) and how these will subsequently evolve into GC's have been described, for example, 
by Elmegreen (1993), Boily, Clarke & Murray (1993), and Murray & Lin (1993). 

HST searches of interacting and merging systems have been spectacularly successful in finding 
compact blue objects which are likely to be YPC's formed during the interactions and mergers. 
Large numbers of blue compact objects were detected in the interacting galaxy NGC 4038/39 (Whit- 
more &; Schweizer 1995), and merger remnants NGC 1275 (Carlson et al. (1998), Faber (1993), 
Holtzman et al. (1992)) and NGC 7252 (Whitmore et al. 1993). The brightest of these clusters 
have the attributes expected of proto-globular clusters (Whitmore 2001). 

Underlying all of this is the hope that, as the merging systems settle to form an elliptical 
galaxy, the YPC's will survive to contribute to its rich globular cluster system (GCS ). An impor- 
tant observational element is missing from this picture: are the kinematics of the YPC's at the 

present time consistent with their ultimate spatial distribution in the cluster system of the final 
elliptical remnant? To answer this question, we need to compare velocity data for the YPC's (or 
for their progenitors) with the kinematics of associated features in the hydrodynamical simulations 
of merging spirals. 

To make this comparison, GC progenitors need to be identified. Assuming YPC's evolve 
into GC then some possible sites for GC formation during mergers include giant H II regions 
and the detached H I fragments and gas clumps within the extended tidal features found around 
merging galaxies (e.g. Due & Mirabel (1994)). We would like to eventually test whether the spatial 
and velocity distributions of these GC progenitor candidates are consistent with the statistical 
distribution and kinematics of small-scale condensations in the hydrodynamical models of merging 
galaxies. We have collected data relevant to measuring the masses and the spatial and velocity 
distributions of H I and H II regions within galaxy systems that sample the range of galaxy-galaxy 
interaction stages leading to a merger remnant. This paper, and a companion paper by English et 
al. (2002: Paper II), describe observations of one merging system, NGC 3256. 

1.2. Globular Cluster Formation Scenarios 

Searle &; Zinn (1978) proposed that GC's form in small fragments which are subsequently 
accreted by the parent galaxy. Their fragments are presumably the low-mass end of the hierarchy 
of fluctuations that aggregate to form a spiral system. In major mergers, another source of fragments 
are those torn from interacting galaxies; falling back into the potential well of the remnant galaxy. 
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they experience star- formation due to tidally induced shocks, and contribute their resultant GC's 
to the emerging elliptical. This globular cluster formation scenario motivated the neutral hydrogen 
observations of NGC 3256 which will be presented in Paper II. 

The ionised hydrogen observations presented here were motivated by Kennicutt and Chu's 
(1988) suggestion that giant H II regions arc the birthplaces of YPC's like the young globular 
clusters seen in the Magellanic Clouds. We consider the possibility that giant H II region formation 
is enhanced in the disks of interacting spirals during the interaction. The associated YPC's would 
be subsequently redistributed, along with the other components of the two disks involved, as the 
merger of the parent galaxies proceeds (e.g. Schweizer (1995)). This is supported by the clustering 
of YPC's in NGC 4038/39. The youngest star clusters spatially cluster together more than the 
intermediate age clusters (which are apparently associated with tidal tails) and than the old GC's 
(distributed mostly throughout the disk) , suggesting mixing has occurred over several orbits (Zhang 
et al. 2001). This is also consistent with the observation that the clusters in giant ellipticals have a 
roughly similar spatial distribution to the diffuse light distribution. [Although in some ellipticals the 
GC distribution is a little flatter than the light distribution (Harris 1993), in others hke NGC 1399 
the cluster and light distributions are almost identical in slope (Harris k. Hanes 1987).] 

The mass spectrum of GC's in ellipticals results from the initial mass spectrum during the 
cluster formation epoch and the effects of subsequent dynamical evolution. Hence the observation 
that, in elliptical galaxies, the mass spectrum of GC's with masses greater than IO^A^q is fit by a 
near universal power law (Harris & Pudritz 1994) which has the same slope as the giant H II region 
mass spectrum (e.g. Kennicutt et al. (1989); McCall et al. (1996)) is interesting in the context of 
this globular cluster formation scenario. 

1.3. The Specific Frequency of Merger Remnants 

The specific frequency Sn of globular clusters is the number of GC's per unit absolute magni- 
tude My = — 15 of galaxy luminosity (Harris &; van den Bergh 1981). The range of Sn observed 
for elliptical galaxies covers an order of magnitude (Harris 1993) and its value for spiral galaxies 
can be an order of magnitude lower than for cD ellipticals which are at the high extreme. 

Schweizer (1988) proposed that the Sn of a remnant elliptical would be larger than the value 
derived from the sum of the GC's populations of the parent galaxies due to GC's formation during 
mergers. Ashman and Zepf (1993) determined the number of GC's per unit stellar mass, using a 
characteristic M/Ly for each morphological type. They find that normal [i.e. non-cD) elliptical 
galaxies have more than twice as many GC's per unit stellar mass as spirals. Therefore even the 
construction of normal ellipticals via the merger of spirals would require the formation of GC's in 
the merger process. 

Of course this would also require GC's to form in preference to non-cluster stars in the merger. 
While this is theoretically plausible (e.g. Elmegreen & Efremov (1997), and vast numbers of YPC's 



-4- 



axe clearly formed in mergers like NGC 4038/9, there is no direct evidence that this is so. Still 
NGC 1275 provides an example of a system with a high Sn (~27) and Carlson et al. (1998) find 
that it will remain larger than that of the old GC's population even if the majority of the low-mass 
clusters are destroyed. 

We should comment here on an important paper by Forbes ct al. (1997) who examined the 
color distribution of GC's in ellipticals. In galaxies with a high specific frequency of GC's (mostly 
cD ellipticals), they found that the excess GC's are mostly metal-poor. This is not consistent 
with forming the high specific frequency of GC's through mergers of spiral galaxies with near-solar 
metallicity. But it does not argue against formation of YPC's in mergers: this is clearly observed 
to happen. The Forbes et al. result suggest that the bulk of GC formation in cD ellipticals may 
have occurred during the early phase of hierarchical merging, while the mean metallicity was still 
low. That is, the GC's may have formed through the early merging of a large number of relatively 
unevolved fragments, rather than through the later merging of more chemically evolved galaxies. 

1.4. NGC 3256 

NGC 3256 is a nearby (2820 km s~^) merging system which is experiencing a spatially extended 
and highly luminous starburst (~ 3 x 10^^ L0; Sargent et al. (1989)). The 10 /xm luminosity of 
the visible Ha nucleus (~ 2 x 10^° L©) rivals that of Seyfcrt galaxies (Graham et al. 1984). The 
system has another "nucleus" 5 arcsec to the south, obscured by dust but detected in non-optical 
wavelengths such as: 3 cm (Norris &: Forbes 1995); near- infrared K-band (e.g. P. MacGregor (see 
Paper II), (Lipari et al. 2000)); and X-ray (Lira et al. 2002). This system also has a star-forming 
complex in its 'disk', an arc of giant H II regions, and 2 extended tidal tails. The 2 X-ray nuclei, 
along with the diffuse X-ray emission, support the scenario that NGC 3256 is powered by a starburst 
rather than an AGN (Lira et al. 2002). Therefore we tentatively assume that NGC 3256 has two 
cores (one from each parent galaxy) that are on the verge of merging. 

Lipari et al. (2002) argue that NGC 3256 is the result of a multiple merger with the disk star 
forming complex as the third nucleus. Since there are only 2 tidal tails, they suggest that either an 
on-going merger system encountered a third galaxy or that 2 gas rich spirals plus a minor galaxy 
merged simultaneously. These cases, for our order of magnitude estimates, can be replaced by a 2 
galaxy prograde interaction scenario. 

In Paper II we estimate the dynamical timescale associated with the neutral hydrogen compo- 
nent of the tails and compare our observations with numerical simulations of starbursts in galaxy 
mergers which generate elliptical remnants (e.g. Mihos k, Hernquist (1994)). In these models the 
timescale since the last pericentre approach of the parent galaxies through to coalescence is about 
500 Myr, comparable to our estimate of the dynamical time (~ 500 h~^ Myr) since pericentre in 
the observed NGC 3256 system. If the model parent galaxies contain dense bulges, a strong star- 
burst is also produced at about 500 Myr post-pericentre. Although the K-band luminosity profile 
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of NGC 3256 is not yet the r^l^ distribution of a fully relaxed system, the profile is consistent with 
a merger phase in which the cores of two galaxies are about to coalesce (Moorwood & Oliva 1994). 
Subsequent violent relaxation would cause this system to have the structure of an elliptical galaxy 
in its final merger state (Schweizer 1986). 

Hernquist and Boltc (1993) examine the GCS's produced in simulations of mergers of two 
disk galaxies. The surface brightness profile of the model GCS's are also well fit by the observed 
rV4 law. They argue that one of the sites of new cluster formation will be in the inner regions 
of the merger remnant, because a large fraction of disk gas is driven to each disk's centre and 
subsequently the cores of the disks merge (Barnes & Hernquist (1992), Mihos h Hernquist (1994)). 
Since NGC 3256 is approaching the merger stage (see the merger sequence by Toomre (1977)), 
it is thus not surprising that H II region formation is enhanced in the galaxy's disk. Again, the 
inner disk gas would provide the high pressure environment to produce bound clusters according 
to Elmegreen k, Efremov (1997). And, though all may not be bound, Zepf et al. (1999) find several 
hundred compact, blue clusters likely to be young clusters within the 7 square arcsec central region. 

Our photometry (§ 2) shows that, although NGC 3256 has relatively few H II regions compared 
to some other interacting systems, these regions are comparable in flux to about 85 30-Doradus-like 
H II regions (§ 3). Since the 30 Doradus nebula contains the YPC NGC 2070 (Meylan (1993)), 
we label as "30-Dor GHR's" any giant H II regions with luminosities greater than or equal to 
1.5 X 10^*^ erg scc^^ (Kcnnicutt & Chu 1988). We present 2-D spectra (§ 2) used to compare the 
velocity distribution of H II regions in this merging system with numerical models; these extend 
spatially beyond the spectra of (Lipari et al. 2000) (§ 3). We estimate the number of YPC's which 
could be born between the last closest approach of the parent galaxies and the currently observed 
epoch in § 3. We summarize in § 4 that NGC 3256 is an emerging elliptical galaxy that is producing 
a population of YPC's consistent with the specific frequencies of GC's that would be expected for 
elliptical remnants formed via mergers of spiral galaxies (Ashman & Zepf 1993). 

2. Observations, Reductions, and Data Analysis 
2.1. CCD Images 

2.1.1. Observations and Reductions 

On February 25, 1991 we obtained 2 x 300 sec broadband I and 2 x 600 sec narrowband Ha 
exposures of NGC 3256, and E region standard field # 5 (Graham 1982), with the 1 metre telescope 
at Siding Spring Observatories (SSO). The EEV CCD has 832 x 1152 pixels, pixel scale of 0.57 
arcsec, readout noise 4 e" and gain 1 e^ per adu. The seeing was better than 2 arcsec. The central 
wavelength of the narrowband filter is 6594 A, its FWHM 52 A and its equivalent width 24 A. 
At the redshift of NGC 3256, the [NII]A6583A emission line lies in the wing of the filter passband 
where its contribution to the image is negligible. 
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We first corrected the non-linearity of this CCD (< 1.5% at 18,000 adu to 6% at the extrap- 
olated value of 70,000 adu). Subtracting a fit to the overscan region of every image reduced the 
low-level time-dependent horizontal bias patterns caused by electronic pickup. The dark current 
of this CCD is negligible so only bias images were subtracted. The target images were flatfielded 
using dome and twilight flatfields for each individual filter. The processed images of the galaxy were 
shifted to a common position and combined for each filter with an algorithm that rejects cosmic 
rays. Image editing corrected residual bad columns and pixels. 

For the continuum images, we used the I-band, thus avoiding contamination of the continuum 
images by Ha and [Nil] emission lines which are included in the R-band. To scale the Ha and 
I-band images, we used aperture photometry of the target galaxy itself. Assuming that the outer 
regions (i.e. envelope) of the galaxy disk do not contain a significant amount of Ha emission, 
the I-band image was initially scaled such that the I-band envelope matched the intensity of the 
envelope in the narrow-band Ha image. We then experimented with scaling I to Ha continuum 
within the inner 11 arcsec radius. We adopted the scaling which avoided negative Ha flux in this 
H II region zone when the scaled I-band image was subtracted from the narrow-band image to form 
a continuum-free Ha image. This is likely to underestimate the true Ha flux. 

The slit spectra of the H II regions have an equivalent width in Ha > 80 AA confirming that 
they are bright. The spectra also indicate that our underestimate in Ha photometry is < 15 %. 

2.1.2. Flux Calibration 

The Ha flux density of our narrow band images was calibrated (in ergs s~^ cm~^) using 
observations of 5 standard stars in the #5 E-region, with the zeropoints given by (Bessell 1990, 
1992). The absolute flux density through the narrowband filter was determined by interpolating 
the broadband zeropoints to the central wavelength of the filter, and using the equivalent width of 
the filter to calculate the integrated continuum light for each standard star over the bandpass. 

We used a circular aperture of 4 arcsec to measure the observed sky-subtracted count rates 
for each star. These were corrected for atmospheric extinction using the mean airmass of the 
observations. The mean extinction coefficient for the Ha filter, fcesge = 0.085, was estimated 
by interpolation between the mean broadband coefficients determined at SSO by Anja Schroder 
(private communication). The estimated error of this calibration is about 6%. 

2.1.3. H 11 Region Photometry 

Each irregularly-shaped giant H II region complex in NGC 3256 was outlined with a polygon. 
The sky-subtracted count rate, determined from the annular sky aperture marked on Figure 1, was 
corrected within each polygon for atmospheric extinction and for Galactic absorption. The Milky 
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Way absorption (Agsgg = 0.35 mag) was determined using the colour excess for NGC 3256 (Burstein 
&: Heiles 1982) and interpolating through average interstellar extinction values A;\/E(B-V) (Savage 
& Mathis 1979). We did not attempt to correct for extinction within NGC 3256 itself. This Ha flux 
was subsequently corrected for the transmission of the filter at the wavelength (6625 A) associated 
with the redshift of the galaxy (2820 km s~^; see determination in § 3.1). The total Ha energy 
output, presented in Table 1, was derived using Hq = 100 h km Mpc~^. 

We did not attempt to divide up large complexes into apparent subcomponents since our goal is 
to estimate whether the complexes are comparable to 30 Doradus in their Ha flux. The compromises 
in polygon shape due to the proximity of other H H regions provide the largest uncertainty in the 
flux estimate (la ~17%), but this is insignificant for our comparisons with the flux of 30 Doradus. 

2.2. Ha Spectra 

2.2.1. Observations 

In order to determine the velocities of the H H regions, long slit spectra were acquired at SSO 
using the Double Beam Spectrograph (DBS) and Photon Counting Array (PCA) at the Nasmyth 
focus of the 2.3m telescope on Mar. 8 & 9, 1991. The red 1200 1 mm~^ grating gave a dispersion of 

0.4 A (~ 19 km s"^) per pixel, a FWHM resolution of about 1 A , and a total wavelength range of 
about 300 A. The slit width was 1.8 arcsec and the seeing about 2 arcsec. In the spatial direction, 
each pixel corresponds to 0.66 arcsec resulting in a useful slit length of about 330 arcsec. 

To place the slit precisely, we placed one end of the slit on one of 3 stars near the interacting 
system (labeled A, B, and C in Fig. 2). Using this star as a pivot point, the spectrograph was 
rotated to various position angles on the sky. For each position angle a number of H II regions lay 
along the slit and hence produced a long-slit spectrum which contained spectra for a number of 
emission-line objects as well as the continuum of the flducial star (which provides a spatial reference 
point) . 

Exposures were typically 1000 sec. A neon lamp was observed before and after every target 
spectrum. This allowed us to correct for the small flexure in the DBS as it rotated to maintain 
constant position angle on the sky. In order to ensure linearity of the detector in the flatfield 
exposures, we took quartz lamp spectra at rates less than 0.1 Hz pixel"^. 

2.2.2. Ha Long- Slit Reductions 

PCA data have no bias, and cosmic rays are rejected in hardware. Therefore, before wavelength 
calibration, it was only necessary to flatfield the images as usual, using quartz lamp spectra: sky 
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fiats provide the correction for vignetting along the sht. 

The wavelength calibrations, determining dispersion solutions and rebinning to linear wave- 
lengths, used standard procedures (Massey et al. 1992). We did not correct for distortions in the 
spatial direction since these were insignificant at the accuracy required by our study. We subtracted 
the background only when night sky OH lines lay very close to the redshifted H II emission lines. 

To measure the radial velocities, wc summed over 3 rows (~ 2 arcec, to match the seeing) 
and determined the centre of the Ha emission-line with a single-gaussian fitting routine. We then 
stepped up one row, again summed 3 rows, and determined the central velocity of the Ha emission 
at this new position. This was repeated for all sections in the 2-D spectrum containing emission. 

The PCA has an intrinsic semi-periodic fringing pattern which changes slowly with time. 
Flatfielding removed most of this pattern. We unsuccessfully attempted to remove the remaining 
fringing by 2D Fourier filtering, and the residual coherent noise ultimately limited the attainable 
S/N and dominated the error in our velocity measurements which increased with a decrease in 
emission intensity. To estimate this error empirically, we produced a calibration curve using the 
sky lines in a set of night sky spectra. A spline curve was fit to a plot of (error in measured 
wavelength) versus (median intensity of sky line). This curve indicated that spectral features with 
peak intensities above 40 counts per wavelength pixel have an rms error of 7 km s^^ (dominated 
by the wavelength calibration error). The error increases with decreasing counts to 18 km s~^ for 
peak intensities of about 10 counts per pixel. 

2.3. Velocity Field Analysis 

An optical velocity field image of H II regions. Fig. 2, was constructed by editing the Ha CCD 
image such that the heliocentric velocity values were substituted at the image pixel coordinates 
associated with the positions of the DBS slit. Gaps in the resultant DBS slit map, due to the 
difference in scale (arcsec/pixel) between the DBS and CCD images, were assigned interpolated 
velocity values. A typical positional uncertainty in the velocity field image is about 1 CCD pixel 
(0.57 arcsec). Due to the long exposure times, the DBS spectra were more sensitive than the 
imaging to diffuse emission. This results in a velocity field data in regions without H II region 
peaks in the image. 

The velocity field image was colour edited to heighten the contrast of consecutive velocity value 
bins. The range of colour bins is nearly linear with each bin spanning 33.5 km s~^. The full colour 

range was chosen to maintain the convention 'red corresponds to redshift and blue corresponds 
to blueshift'.^ To create the intervening velocity range, the red and blue extrema were blended 



^The saturation of the red was set such that it was dull and the blue chosen was a 'warm', bright blue. Hence the 
red colour, in terms of human visual perception, appears to recede and the blue colour visually approaches, in the 
same sense as redshift and blueshift. 



-9- 



together such that the neutral gray bin in Fig. 2 contains the hehocentric systemic velocity value. 

The data from the horizontal slit position (associated with fiducial star A) is also presented 
in plot format (Fig. 3) in order to provide an indication of the galaxy's velocity curve. Since 
the inclination of this galaxy is close to face-on wc have not attempted to deproject the data. 
Even when multiple gaussian features were present in the 1-D spectra, they have been assigned 
an average velocity; the errors given are flux dependent errors (see § 2.2.2) and do not reflect the 
range in velocity that these features span. 

Also in Table 1 we present the heliocentric velocities corresponding to the point on the DBS slit 
which is closest to the centre of each polygon used in the photometry (see Fig. 1). The difference 
between a polygon's centre and the position of the peak intensity of the H II region can be up 
to a few arcsec. Also, the slit positions generally did not pass directly through the centre of the 
H II region. A comparison of the two slits which pass near the northern nucleus indicates that the 
velocities tabulated in Table 1 may differ from the mean velocity of a given H II complex by up to 
20 km s-^ 



2.4. Relationship to Other Data 

These Ha images and spectra have in some respects been superseded by data acquired by 
Lipari et al. (2000). Their Ha image from the ESO New Technology Telescope is also used by Lira 
et al. (2002), who calibrated the flux (to within a factor of 2) of the Lipari et al. image using the 
Lipari et al. long slit spectroscopy of several H II knotts. However Lira et al. only report a total 
flux (~ 2 X 10~^^ erg s~^ cm~^ in an aperture with a 41 arcsec diameter). Lipari et al. (2000) 
measured the Ha flux of H II knotts via spectroscopy, effectively subdividing our selection of H II 
regions without measuring all the flux per region. For example, they present the flux of regions 
labeled Rl, R2, and R4 but not of R8 and R9 although all of these would be in our regions Bl and 
B2; see Table 3. Although our measurements are therefore difficult to compare, their R7 region's 
Ha flux differs from our measurement of B6 by 8%. 

Our spectra are included not only for completeness sake but also because the resolution is 34 
km s~^(with errors ranging from 7 to 18 km s~^) compared to the Lipari et al. (2000) instrument 
resolution of 90 km s~-^ (errors ranging from 15 to 30 km s~^). Additionally Lipari et al. only 
present velocity field maps with a width of 60 arcsec. Researchers, particularly those simulating 
the interactions between galaxies, may find our 330 arcsec long slit provides useful information on 
almost all radii associated with the ionized hydrogen emission component of this galaxy and that 
our velocity field reveals interesting anomalies beyond the inner region (e.g. redshifted emission in 
the tail to the SW in Fig 2.) 
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3. Results 

3.1. Velocity Curve and the Dynamical Mass 

Visual inspection of Fig. 2 shows largcscale rotation, in the same sense as the rotation observed 
in the neutral hydrogen data (Paper II). However the velocity behaviour of the optical disk appears 
less disturbed by the merger activity than do the H I arms. 

The peculiar morphology of NGC 3256, the uncertainty of its centre, the position angle of the 
major axis, and its near face-on inclination, make it difficult to extract an accurate deprojected 
circular velocity curve. Fig. 3 shows a velocity curve from a spectrum with the slit at position angle 
(P.A.) 90° . Examination of Ha features in Fig. 2 shows that the turnover points in the velocity 
curve are representative of the velocity range of the system. The value at the midway point between 
the turnover peaks is 2820 it 11 km s~^ and the velocity difference, AV, gives a rough estimate of 
the projected rotation velocity; Vdrc = AV/2 107 it 11 km s~^. 

This agrees with theFeast and Robertson (1978) estimate of the systemic velocity and rotation 
amplitude (using P.A. = 100°), Lipari et al. (2000) (using P.A.= 90°, however only for an inner 
region 40 arcsec square), and with the H I position- velocity (PV) diagram shown in Paper II. 
(We note that the distribution of H I in the PV diagram means that the rotation amplitude is 
significantly larger than the separation of the H I horns in the integrated H I profile also shown in 
Paper II.) 

Assuming that the rotation curve is fiat, we estimate that the mass within radius r is M{r) = 
3.2 X 10^ [vcirc(km s"-*^)]^ r(arcsec) h^^ Mq. This is no more than an order of magnitude estimate 
because the system is unlikely to be in centrifugal equilibrium at this stage. Dynamical mass limits 
are given in Table 2 for a few radii associated with structural or kinematic features. 

3.2. H II Regions and Globular Clusters 

Our goal in this section is to consider whether the H II regions observed in NGC 3256 could be 
GC progenitors. For an ionised hydrogen region to be a GC progenitor, in the sense of Kennicutt 
and Chu (1988), we expect it to have an Ha luminosity of at least 1.5 x lO^'^ ergs sec~^. This is 
comparable to the total Ha luminosity of 30 Doradus and corresponds to an ionised hydrogen mass 
of about 8.5 x 10^ Mq (Kennicutt 1984). We refer to these giant H II regions as 30-Dor GHR's. 
Additionally, each 30-Dor GHR's should contain a rich star cluster; 30 Doradus harbours R136 
which, from HST images (e.g. Meylan (1993)), is known to have the the dense core and extended 
halo structure typical of a young globular cluster or YPC. Since the timescale for the formation of 
a YPC is comparable to the lifetime of a giant H II region (Kennicutt &: Chu 1988), it is highly 
likely that the cores of objects like 30 Doradus will be identified as YPC's when the surrounding 
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OB associations have faded. 

HST observations of NGC 4038/9 (e.g. Whitmore and Schweizer (1994)) show both spectacular 
populations of giant H II regions and hundreds of YPC's . Evidence of an association between 
YPC's and 30-Dor GHR's is provided by a comparison of YPC's with multiwavelength emission 
in NGC 4038/9 by Zhang et al. (2001). They find that red clusters (age < 5 Myr) are associated 
with radio continuum, CO and Infrared emission, that blue young clusters (3 < age < 16 Myr) 
are correlated with Ha and Xray emission, and bright, older clusters (16 < age < 160 Myr) are 
correlated with Xray emission. They note that this is consistent with the cluster heating the dusty 
environment it formed in, then ionizing the surrounding gas, and finally expelling the surrounding 
material via stellar winds and super novae. 

In NGC 3256 Lira et al. (2002) observed diffuse Xray emission, using Chandra Xray Obser- 
vatory, that spatially co-incides with diffuse Ha emission and they claim that this is evidence that 
the Xray gas is heated by a starburst (rather than an AGN). Additionally most of their discrete 
Xray sources are co-incident with H II knotts; see Table 3 for correlations with our observations 
of diffuse emission. The potential sources generating the discrete Xray emission include supernova 
remnants and Xray binaries. This suggests the sources are in clusters with massive stars, which 
may be initially embedded in H II regions in a manner consistent with the Zhang et al. picture. 

This picture is echoed by Alonso-Herrero et al. (2002) who specifically look for a spatial corre- 
spondence between H II regions and YPC in NGC 3256 using NICMOS Paa and Infrared images. 
They use evolutionary synthesis models to indicate that the ~8% coincidence rate, compared to the 
number of H II regions + star clusters, is due to evolution effects and their detection thresholds. 
That is, after 9 Myr the H-band luminosity, used to detect YPC, peaks while the Ha luminos- 
ity drops below their detection limit. Since the H-band luminosity threshold allows detection of 
clusters from ages 0-100 Myr, more YPC are detected than H II regions. Those H II regions with- 
out detected IR clusters probably harbour clusters which are obscured by dust in the manner of 
Galactic H II regions less than 5 Myr old. Alonso-Herrero et al. (2002) believe that they are only 
detecting spatial coincidences when the YPC have masses ~ W^Mq and ages less than 7 Myr. 

At optical wavelengths Zepf et al. (1999) find about 1000 blue, compact clusters in the inner 
region of the NGC 3256 merger. Therefore, at an order of magnitude level, we can compare the 
number of clusters like R136 to an estimate of the specific frequency of GC's required if NGC 3256 
is a merger of 2 galaxies. 

To determine the number of R136-like clusters we first select only those Zepf et al. clusters 
with -0.55 < (B-I)< 0.5. Bruzual and Chariot (1998) models generate unique ages between 3 and 
7 Myr for clusters in this color range. As well, models by Schaerer and de Koter (1997 and papers 
listed therein) suggest O and B stars provide ionizing radiation for about 7 Myr, making these 
clusters comparable in age to the observed H II regions which might harbour them. 
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The observed absolute magnitudes (Mb) of these blue YPC's are calculated using Hq = 75 
km Mpc^-^. The change in the Mb of R136 (-12.6; Kennicutt & Chu (1988)) is also calculated 
for the 3 to 7 Myr age range according to the Bruzual and Chariot (1998) models. Examining the 
difference between the observed magnitude of a cluster and the magnitude predicted for R136, at 
same age as the cluster, indicates that 17 out of 196 of the bluest clusters have luminosities, and 
hence masses, greater than R136. The total mass in the 17 clusters is equivalent to about 45 x the 
mass in R136 (AIrisg)- (The Alonso-Hcrrcro ct al. (2002) clusters that arc coincident with H II 
regions probably have ages between 3-7 Myr and masses W^Mq. At least 18 coincidences occur 
in the 19.5 square arcsec field-of-view of NIC2. Scaling this area to match the 36 square arcscc area 
of the PC chip predicts the existance of 60 IR clusters with R136-like masses, which is comparable 
to our estimate of 45.) Extrapolating to the ~1000 YPC's detected by Zepf et al. (1999) suggests 
that there are roughly 100 YPC's with masses greater than R136 or a total of 225 Mrisg clusters. 

Thus we estimate the existence of 100 - 200 massive YPC's; 100 YPC's if the masses are 
distributed amongst clusters in a manner similar to this sample and approaching 200 YPC's if 
more of the clusters are about 1 A1ri36- This is a conservative estimate since only a limited field of 
view has been studied. Additionally our threshold of (B-I)< 0.5 eliminates blue clusters reddened 
by intervening dust and we note that all 17 clusters are situated in the less dusty east section of 
the disk. 

An alternative YPC estimate uses Ha photometry and also suggests that more than 200 R136- 
like clusters could form over the lifetime of the galaxy-galaxy interaction. 5 H II complexes in 
Table 1 have fluxes above 1.5 x 10^° ergs sec~^, giving a total of 85 30-Dor GHR's not including 
those that might exist in the obscured southern nucleus. (This is consistant with Alonso-Herrero 
et al. (2002) who detect at least 50 30-Dor GHR's in the FOV of NIC2, which includes only Bl, 
B5, part of B6 and B8. Our remaining H II complexes contribute 67% of the Ha luminosity due to 
30-Dor GHR's. This implies about 35 30-Dor GHR's were not imaged by NIC2.) All the Zepf et al. 
clusters but one project onto the sky in the vicinity of the 3 30-Dor GHR's regions in the 7 kpc by 
7 kpc field of view that they analyzed. The number of bright, blue R136-like clusters associated 
with the 30-Dor GHR's suggests that 20-50% of these 30-Dor GHR's may currently cradle YPC's. 
Notably, Alonso-Herrero et al. (2002) claim about 55% of their high-luminosity H II regions have 
an IR cluster counterpart. 

To estimate the total number of YPC's which might form from 30-Dor GHR's over the era of 
the current interaction stage, we consider the rate of 30-Dor GHR's formation in NGC 3256. Using 
a lifetime of 7 Myr gives a current rate of 12 Myr~^ (compared to less than 0.2 Myr~-^ for a typical 
Sc galaxy with Mb = -19.5 Kennicutt & Chu (1988)). This enhanced rate could be maintained 
for longer than the age of the nuclear starburst (between 10 and 27 Myr; Doyon et al. (1994) and 
Rigopoulou et al. (1996)). However the numerical merger simulations of Mihos &; Hernquist (1994), 
in which each parent galaxy has a bulge, suggest that this rate has probably not been continuous. 
The global star formation rate of each parent galaxy increases by only a few percent until the 
galaxy cores coalesce, and then the star formation rate increases dramatically. 
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Assuming that the cores are now about to merge, the time from pericentre through to coales- 
cence is about 50 model units, corresponding to about 500 h~^ Myr in NGC 3256 (see Paper II). The 
peak of the star formation extends over about 4 model units, or about 40 h"-*^ Myr in NGC 3256. 
If we are now observing the peak star formation rate, then this time interval generates at least 
640 30-Dor GHR's (for h = 0.75). The formation rate for the previous 400 h~^ Myr is about 
7% of the peak rate value. Hence another 250 H II regions probably formed between the time of 
pericentre separation and the observed enhanced star formation epoch. If 20-50% of the 30-Dor 
GHR's harbour R136-like clusters, then we might expect at least 200-400 YPC's to form over the 
lifetime of the interaction. 

Although there is some controversy about the relationship between the YPC mass function 
and the GC mass function (see Fritze - v. Alvensleben (2001) and papers therein), we believe that 
the R136-like clusters in the Zepf et al. data will acquire GC colors and luminosities over several 
Gyr. It is difficult to estimate how many of the R136-like clusters would survive the evolution 
of the NGC 3256 GCS over the next several Gyr. However if the current power-law luminosity 
function of the YPC system (Zepf et al. 1999) transforms to the log-normal distribution of an 
older GCS mainly via evaporation (due to internal relaxation) then the loss of mass occurs mainly 
among clusters less massive than those in the median mass bin (Fall &: Zhang, in press). Since the 
mass of R136 corresponds to the median bin values, we expect that these couple of hundred YPC's 
may survive. 

If typical ellipticals form via the merging of spirals, the number of new massive clusters formed 
in the merger needs to be at least comparable to the number of globular clusters originally associated 
with the spirals {e.g. Ashman & Zepf (1993)). Hence we now estimate the original number of 
GC's and compare it with the number of YPC's in NGC 3256. To get an order of magnitude 

estimate of the number of GC's associated with NGC 3256's parent systems, we arbitrarily assume 
that the parent galaxies have similar bulge/disk ratios to M31, and compare the dynamical mass 
(~ 3 X lO^'^ h^^ -^0) within the radius of the I-band envelope of NGC 3256 with the dynamical 
mass of M 31 (3 X 10-*^^ M.q) within the 20 kpc radius of its globular cluster system (GCS) (Huchra 
1993). Since M 31 has about 400 confirmed GC's (Fusi Pecci et al. 1993), scaling by the ratio of 
dynamical masses suggests that the original number of GC's belonging to the two-parent system 
was roughly 40-80. The larger value uses the Feast and Robertson (1978) inclination and position 
angle (see §3.1). As a check we use the specific frequency per unit mass for disk galaxies (2.2 
clusters x 10^ Mq^; Zepf & Ashman (1993)) which gives 66 GC's . 

In summary, since ellipticals have typically more than twice as many GC's per unit mass as 
spirals (e.g. Ashman & Zepf (1993)), we need at least as many GC's to form in the merger as 
originally existed in the parent system. If NGC 3256 is the result of the interaction of 2 similar disk 
galaxies, as indicated by its 2 optical tails, then the number of GC's required would be less than 
100. Extrapolating from the fraction of the bluest YPC's observed in NGC 3256, roughly 100-200 
YPC's are expected to have masses equal to or greater than R136. (Using Ha photometry and 
numerical models suggests 200-400 R136-like clusters form.) These should fade to GC luminosities 
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over several Gyr and many are expected to survive destruction, since it is the lower mass clusters 
which evaporate. Hence the observed YPC's appear numerous enough to populate the GCS of the 
emerging galaxy such that its specific frequency will be consistent with that observed in typical 
elliptical galaxies. 

Additionally, the bluest, massive YPC's are located in the vicinity of observed 30-Dor GHR'ss, 
contributing to the notion that some fraction of 30-Dor GHR's do cradle massive YPC's , as 30 
Dor harbors R136. If interactions induce the formation of 30-Dor GHR's, the observed luminosities 
indicate that almost 900 30-Dor GHR's would form in NGC 3256 throughout its merger epoch. 
In order for 30-Dor GHR's to be considered GC progenitors, the specific frequency argument 
only requires that about 10% of these harbor YPC's which survive several Gyr and have masses 

> A^R136- 

4. Conclusions 

An approach, outlined in our introduction, that would be fruitful for assessing whether giant 
H II regions may be globular cluster birth places, uses hydrodynamic modeling in order to see 
if the kinematics of the H II regions at various interaction and merger stages is consistent with 
the final spatial distribution of GC's in elliptical (i.e. merger remnant) galaxies. This approach 
should be applied to a sample of galaxies spanning the interaction stages that lead to a merger 
remnant. It requires a comparison of observed H II region positions and velocities throughout the 
merger sequence with the statistical position and velocity distributions of condensations of stellar 
and gas particles seen in simulations which generate elliptical-like merger remnants (e.g. Mihos &: 
Hernquist (1996)). At every interaction stage, the behaviour of the observed progenitor candidate 
should be similar to that of the condensations in the model. If this picture is correct both the 
observed and numerical sequence should produce a GCS that is consistent with GCS's observed 
around elliptical galaxies. 

Another approach compares the observed specific frequency of GC's in typical ellipticals with 

an estimate of the number of GC's formed as two similar disk galaxies interact and coalesce to form 
an elliptical merger remnant. Since ellipticals have typically twice as many GC's per unit mass 
as spirals, at least as many GC's need to be created in the merger process as originally existed in 
the sum of the 2 parent galaxies (Ashman &; Zepf 1993). The argument is that tidal disturbances 
provide a mechanism for enhancing star formation and hence generate giant H II regions in the disk 
of the emerging elliptical. On theoretical grounds interactions are expected to lead to enhanced 
H II region formation (e.g. Larson (1987)), and observations demonstrate that giant H II region 
populations are enhanced at various stages of the interaction-through-to-merger sequence (e.g. 
NGC 4038/9, NGC 5426). These ionised complexes in turn "cradle" YPC's. (Support that YPC's 
are associated with H II regions is provided by their correlation in NGC 4038/39 (Zhang et al. 
2001).) Massive YPC's are expected to survive evaporation from the GC system of the emerging 
elliptical and we believe that these YPC's will fade over several Gyr such that they would be 
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recognized as GC's.The optical data on NGC 3256 presented in this paper allows us to explore this 
approach. 

To get an order of magnitude estimate of the number of GC's associated with NGC 3256's 

parent systems, we arbitrarily assumed that the parent galaxies have similar bulge/disk ratios to 
M31. Using M31's number of GC per unit mass suggests the original number of GC's belonging to 
the two-parent system was roughly 40-80. Thus over NGC 3256's merger period about 100 YPC's 
must form and survive, if the above scenario is correct. 

Although NGC 3256 has relatively few (7) H II complexes compared to some other interacting 
systems, 5 of these regions are comparable in flux to about 85 30-Dor giant H II regions (30-Dor 
GHR's; each with an Ha flux of 1.5 xlO^° erg sec"^). Given that these 30-Dor GHR's formed 
within the previous 7 Myr, this merging system currently has a giant H II region formation rate 60 
times that observed in an Sc galaxy with Mb = —19.5 (Kcnnicutt & Chu 1988). The 30 Dor nebula 
contains the YPC R136 within NGC 2070 (Meylan 1993) and the formation timcscalc for YPC's 
is comparable to the lifetime of a giant H II region (Kennicutt k. Chu 1988). As well, clusters as 
massive as R136 are expected to survive evaporation due to internal relaxation (S. M. Fall, private 
communication) and hence become part of the remnant elliptical's GCS. Therefore we used this 
flux as a minimum criterion for identifying which H II regions are potential birth places of YPC's 
(Kennicutt & Chu 1988). 

An order of magnitude estimate of the total number of YPC's created throughout the merger 
sequence experienced by NGC 3256's parent galaxies was determined 2 ways. For one estimate we 
determined that about 1/10 of the '~ 200 bluest clusters observed in the inner region of NGC 3256 
by Zepf et al. (1999) have luminosities, and hence masses, greater than R136. Extrapolating 
to the ~ 1000 detected clusters implies 100-200 YPC's exist. The alternative estimate used Ha 
photometry and the variation in star formation rate over the interaction-through-merger timescale, 
calculated using existing numerical simulations (Mihos & Hernquist 1994). These indicated about 
900 30-Dor GHR's formed since the last pericentre approach of the parent galaxies, yielding 200- 
400 YPC's if 20-50% of 30-Dor GHR's have R136-like clusters. We expect clusters this massive to 
survive destruction. 

If NGC 3256 is typical of the merging systems that form globular clusters, and if about 10% 
of the 30-Dor GHR's harbour clusters that survive several Gyr, this would sufficient to explain the 
high specific frequencies of clusters in ellipticals. 

We are grateful to Mike Dopita for suggestions on observing strategy, to Peter Quinn for per- 
spective on interactions and mergers, to Peter Wood for the use of his CCD image linearizing script, 
and to the technical staff at Mount Stromlo and Siding Spring Observatories. J. E. acknowledges 
the support of an Australian National University Postgraduate Scholarship. 
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Fig. 1. — Photometry of H II regions. The underlying Ha image is continuum subtracted. The 

annulus marks the region from which the mode sky value per pixel was determined. The polygons 
mark the areas measured. Those designated Bl through B7 are treated as H II regions and are 
described in Table 1. Comparisons with positions of Ha emission knotts and X-ray sources in 
other projects are in Table 3. The area of B8 lies close to an optically-obscured radio source. The 
inner radius of the annulus is 43 arcsec. 

Fig. 2. — The velocity field of the ionised gas in NGC 3256. The underlying image is the ionised 
hydrogen data acquired using the SSO 1 metre telescope. The spectroscopic data was acquired using 
the SSO 2.3 metre telescope, the Double Beam Spectrograph (DBS), and the Photon Counting 
Arrays (PC A). Each 'line' consists of the velocity values at a slit position. Red indicates receding 
velocities; blue indicates approaching velocities. The construction of this map is detailed in § 2.3. 
The stars used as fiducial pivot points for the DBS slit positions are labeled. 

Fig. 3. — Ionised hydrogen velocity curve. The DBS slit at position angle 90° and anchored at star 
A; see Figure. 2. The declination of this observation is north of northern nucleus. The rms error 
in velocity, determined using the calibration curve described in § 2.2.2, is 7 km s^^ between lOh 
27m 49.7s and lOh 27m 53s. Outside this R.A. range the counts decrease and the error increases 
to 18 km s"-*^. 



Table 3. Designations of Ha emission measurements plus discrete X-ray sources. 





Ha Emission 




X-Ray Sources 


Figure 1 


Feast & Roberston 


T • • J 1 

Lipari et al. 


T • J 1 

Lira et al. 


Bl 


T 


Rl, R9 


7,6 


B2 


Y, Z 


R2, R8, R4 


10 


B3 


X 


R5 


11 


B4 


w 


R15 


9 


B5 




Rll, RIO, R6 


(3 is just south) 


B6 


s 


R7 


2 


B7 




R19 





Note. — See Figure 1 for H II region designations for our data. 
Note that our measurements do not precisely coincide with the H II 
emission regions selected by Feast & Robertson (1978). Also Lipari 
et al. (2000) measure smaller Ha knotts and do not account for all 
of the emission in our selected areas. We only include the discrete 
X-ray sources plotted in the figures of Lira et al. (2002). Besides 
the southern nucleus, Xray sources 5 and 12 are missing from this 
table; our Ha spectra mapped in Fig. 2 suggest 5 is associated with 
diffuse emission while 12 occurs in an apparent gap in Ha emission. 



